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SUMMARY 

Data are presented on loads and motions of a NASA F-106B airplane flying inside 
thunderstorms; the data cover operations in areas of thunderstorms that produced air- 
borne weather radar return between 29 and 50 dBZ within 200 miles of the Langley 
Research Center. 

Three pilots were involved; no significant differences in piloting techniques 
were observed. The thunderstorm flight data were compared with data from other 
thunderstorm missions and from three airline transport airplanes over their 'total 
passenger-carrying operational spectrum from takeoff to landing. The data compari- 
sons indicate that airliners in normal operations occasionally encounter turbulence 
almost as severe as that encountered in these thunderstorm flights. The maximum 
derived gust velocity U^g of ±45 ft/sec measured in this program was about two 
thirds of the 65 ft/sec measured in a Thunderstorm Project conducted in the early 
to mid-1950's. The lower U<j e experienced may be a result of the restriction in 
thunderstorm intensity imposed in this project. 


INTRODUCTION 

The National Aeronautics and Space Administration has been conducting a study 
of the hazards involved in operating aircraft in severe storms. Data are being 
obtained by flying a highly instrumented F-106B airplane through thunderstorms 
within 200 miles of the Langley Research Center, Hampton, Virginia, as described in 
references 1, 2, and 3. Current emphasis is on lightning hazard characterization, 
and over 600 direct strikes to the airplane have been experienced. In addition, the 
airplane loads, motions, and control inputs experienced in the thunderstorm environ- 
ment have been examined and are given in this report. These data have been compared 
with results from 1941 through 1947 "Thunderstorm Project" (ref. 4) , the "Rough 
Rider" project (ref. 5) , and NACA/NASA "VGH" and "Digital VGH" programs taken during 
normal operations of commercial air transports (refs. 6 and 7) . Although in normal 
transport operations, thunderstorms are avoided when possible, sometimes they are 
unavoidable and in other instances the airplane is inadvertently flown into a 
thunderstorm. Therefore knowledge of that environment is needed. Methods of ana- 
lyzing gust loading of airplanes used herein are presented in references 7 
through 10. 


SYMBOLS 

a n incremental normal acceleration from steady-state lg unaccelerated flight, 

g units 

c mean wing chord, ft 

dBZ logarithmic scale of radar measure of water content 

2 

acceleration due to gravity, 32.2 ft/sec 
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TEST EQUIPMENT 
F-106B Research Airplane 

A highly instrumented F-106B "Delta Dart" airplane (fig. 1) is used to make 
thunderstorm penetrations in the Storm Hazards Program. Details on the F-106B air- 
plane and the criteria used in choosing the airplane for this mission can be found 
in references 1 and 2. Prior to each thunderstorm season the airplane is prepared 
to safely take direct lightning strikes (lightning hardened) and the preparation 
verified by a ground test simulating lightning strikes. The lightning hardening pro- 
cedures and verification tests are described in references 1 and 3. Although the 
airplane has undergone extensive modifications to incorporate the instrumentation 
system, the aerodynamic characteristics remain unchanged. Characteristics of the 
airplane pertinent to the evaluation of the data are as follows: 


Takeoff weight, lb 

Typical fuel burn in 90-min flight, lb 

Wing area, S, ft 2 

Fin and rudder area, ft 2 

Wing mean geometric chord, c, ft .... 
Lift-curve slope, m, per radian 


35 000 
6000 
700 
105 
23.76 
2.9 


Airborne Instrumentation Systems 

The F-106B is equipped with a number of data systems to measure the environ- 
mental and electromagnetic characteristics of thunderstorms during penetrations. The 
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Aircraft Instrumentation System (AIS) measures and records all the airplane environ- 
mental and performance data. Accelerations were sensed by the accelerometers 
mounted as close to the center of gravity as practical such that the measurement 
error due to displacement from the center of gravity will be no greater than 
2 percent. 


Test Procedures and Data Selection 

All the pilots agreed to fly the thunderstorm penetrations at a nominal indi- 
cated airspeed of 300 knots. Usually, ATC assigned a fixed altitude and the airplane 
was maneuvered to adhere to that altitude. Also, some maneuvering was included when 
a line of several cells was penetrated. The cells usually did not line up in a 
straight line and sometimes a turn would be initiated before a cell was exited and 
sometimes the next cell would be entered before the turns were completed. The data 
used for this paper were truncated at the beginning and end of penetrations to mini- 
mize maneuver influences on the data shown. All data shown were taken inside the 
thunderstorms . 

Figure 2 shows the number of penetrations versus penetration duration during 
the summer 1982 season; 192 thunderstorm penetrations were involved with a total time 
inside thunderstorms equal to 398.1 min (6.63 hr). The average penetration lasted 
about 2.07 min. The dates, places, and flight number for each pilot is given in 
table I. Flights analyzed in this report comprise only 16 percent of the total num- 
ber of flights in the program through 1984. The airborne weather radar return for 
the area of the storms that the airplane penetrated was between 29 and 50 dBZ, 
although the maximum radar reflectivity frequently exceeded the 50 dBZ operating 
limit. Storm return on the airborne radar was used for guidance into the storm; 

29 dBZ was the lowest level the radar would show. Penetrations into areas of the 
storm with levels in excess of 50 dBZ were avoided to preclude hail encounters. 

Figure 3 shows the distribution of penetration altitudes which range from 10 000 to 
38 000 ft pressure altitude. 

The level crossings of various flight parameters were counted when the value 
crossed the given level with a positive slope for positive levels, and when it 
crossed the given level with a negative slope for negative levels; zero was con- 
sidered a positive level (ref. 6) . The same level can be counted more than once 
without the value becoming zero. Also, the maximum positive and negative values 
were established for each penetration; each curve dealing with maximum positive or 
negative value therefore represents 384 data points. The derived gust velocity was 
calculated by using equation (13) from reference 8; that is. 


2Wa n 

Ude = K g p o V e mS 


with a constant lift-curve slope of 2.9. Fuel burn was taken into account to 
establish airplane weight; weight decrease due to fuel burn was about 6000 lb during 
a typical flight lasting 90 minutes. During storm penetrations, changes in normal 
acceleration were taken with respect to steady-state lg conditions; altitude changes 
were taken with respect to the average altitude for that penetration. 
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RESULTS AND DISCUSSION 


Data Presentation 


The histogram of penetration duration in figure 2 shows that most penetrations 
had durations from 0.5 to 2.5 min; the longest durations were up to 9 min. Since 
the data from some penetrations were truncated to eliminate maneuvers, the actual 
time in the thunderstorm was somewhat higher than shown in this figure. Although 
not a rule, generally the longer the penetration time the bigger and more intense 
was the thunderstorm. Duration thus gives a qualitative indication of the distribu- 
tion of storm intensity. The data in figure 3 indicate that about 68 percent of the 
total time was spent in the pressure altitude interval from 24 000 to 32 000 ft, 
with the total range being from 10 000 to 38 000 ft. 


In figures 4 through 12, data are presented in four forms for each of the fol- 
lowing parameters: normal acceleration, lateral acceleration, longitudinal stick 

position, pitch attitude, lateral stick position, roll attitude, changes in pressure 
altitude, indicated airspeed, and derived effective gust velocity. In each figure, 
the total level crossing results are presented in part (a) ; the level crossings for 
each of the three pilots are presented in part (b) . The maximum positive and nega- 
tive values of each parameter in each penetration are presented as a percentage 
basis for all data in part (c) and for each pilot in part (d) . 


Normal acceleration data are shown in figure 4. A straight line faired through 
the level crossing rates in figure 4 (a) will intercept the a n = 0 line at about 
3000 counts/hr or N Q « 0.83 upward crossing/sec. The equations needed to calcu- 
late N Q by means of the methods of the theory of reference 9 are 


where 


k Q = f (y) 

was obtained from figure 9 of reference 9 and 


y 


= y r 


2W 

mpcgS 


Using these equations for the average case of 28 000-ft altitude and 32 000-lb air- 
craft weight gives y equal to 39.02, and from figure 9 of reference 9, k Q is 
equal to 0.084. The penetration velocity at an altitude of 28 000 ft was 760 ft/sec 
which will yield N Q = 0.855 upward crossing/sec or 3079 counts/hr which agrees 
well with the approximately 3000 counts/hr indicated in figure 4(a). For the con- 
ditions of these flights, that is, flight inside storms with a radar reflectivity 
of no greater than 50 dBZ, the largest excursion was 1.5g (fig. 4(a)), and for 
about 90 percent of the penetrations the maximum excursions were less than ±lg 
(fig. 4(c)). Lateral acceleration data are given in figure 5 with excursions 
generally being less than ±0.3g. About 90 percent of the penetrations had maximum 
lateral excursions less than ±0.2g (fig. 5(c)). The normal and lateral acceleration 
results shown in figures 4(b) and 5(b) show some differences for load increments for 
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flights flown by the three different pilots. Pilot 1 experienced somewhat higher 
normal and lateral acceleration increments (figs. 4(b) and 5(b)) than pilot 2, with 
usually slightly less control input activity (figs. 6(b) and 8(b)). Pilot 3 had 
only a few flights resulting in a small data sample and considerable scatter in the 
data points shown for him. These observations and the U de results shown in fig- 
ure 12(b) may indicate pilot 1 encountered slightly more severe storms than pilot 2. 

Attitude upsets due to storm turbulence and corresponding control inputs to 
counter the storm effects are indicated by figures 6 through 9. Comparison of data 
in figures 7(a) and 9(a) show that an attitude excursion of 10° in roll is almost 
three orders of magnitude more likely to occur than a 10° pitch attitude excursion 
(relative to the trim pitch attitude of about 4°). The altitude excursions shown in 
figure 10 indicate an equal distribution of updrafts and downdrafts. Figure 10(a) 
shows that excursions of ±1000 ft were encountered; figure 10(c) indicates that alti- 
tude deviations less than 150 ft were found in about one third of the penetrations. 
The indicated airspeed variations (fig. 11) reflect primarily the gust component 
parallel to the direction of flight. These data (fig. 11(b)) show good consistency 
between the three pilots. 

Figure 12 (a) shows derived gust velocities U de up to ±45 ft/sec were experi- 
enced by the F-106B but over 50 percent of the penetrations resulted in values of 
U de less than ±16 ft/sec (fig. 12(c)). The curve shown in figure 12(a) for positive 
U de represents the faired average curve that is obtained by considering the negative 
values jointly with the positive values in an absolute value sense. Figure 12(b) 
shows distributions of U de for each pilot. 

Previous studies (i.e., ref. 9) indicate that an equation of the following type 
usually fits the storm data of figure 12(a) quite well: 


N 


PN e 
o 


-Ude/°d 


where P is the proportion of time spent in turbulence, N Q represents upward zero 
crossings per unit of time, and o d is the root-mean-square (rms) value of the 
turbulence intensity. The data of this report also support this equation. The 
specific equation representing the faired curve in figure 12 (a) is 


N 


3200e 


-U de /4.79 


Since the flights were virtually all in turbulence, P = 1 in this study, the equa- 
tions then indicate that N Q = 3200 counts/hr, or N Q = 0.889 counts/sec, and 
a d = 4.79 fps. This value of c? d is a derived equivalent gust velocity at sea 
level. To establish true airspeed values both U de and O d should be multiplied 


by 



p o is sea-level density, and 


p is the density at altitude. An 


average altitude of 28 000 ft is assumed here; thus. 



1.576 
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The true rms value is thus found to be a ^ rue = 7.55 fps. This agrees well with the 
rms values found in previous flight studies through cumulus clouds only, 7.82 ft/sec 
from figure 2 of reference 10, but is significantly less than the maximum 0 true of 
16 ft/sec for flight in severe storms in reference 10 (p. 9) . The restriction to 
dBZ levels below 50 in the present program undoubtedly accounts for the lack of U^ 
approaching the 16 ft/sec measured in the earlier program. 


Gust Data Comparisons 

Other thunderstorm project results .- In figure 13, the derived gust velocities 
for the F-106B are compared with the data from reference 4 for the Thunderstorm 
Project and the thunderstorm data from reference 5 for the Rough Rider project. 

These data show upward gusts only. The data for the F-106B show general agreement 
with the previous storm data up to U^g = 45 ft/sec, the maximum measured, but the 
Thunderstorm Project described in reference 4 measured values of U de approaching 
65 ft/sec. The Rough Rider data shown in figure 13 have been extrapolated in refer- 
ence 5 to higher derived gust velocities than were actually measured. 

Air transport results .- The latest a n and U<j e data from three different air 
transport airplanes were available from the NASA Digital VGH (DVGH) program and were 
in the same format, from the same data reduction methods, so that direct comparisons 
were immediately available. The complete DVGH data set is at present unpublished, 
although the program is described in references 6 and 7. Over 1600 hours of data 
were collected for each of the three types of airplanes: Boeing 727, Boeing 747, 

and Lockheed L-1011. The 727 and 747 flights were in the continental United States 
and southern Canada with some 747 flights being transatlantic. The L-1011 flights 
were in the eastern United States and Puerto Rico. Data were taken from 15 sec after 
takeoff until 15 sec before touchdown. The data were numerically filtered to sepa- 
rate maneuver-induced accelerations and gust-induced accelerations. Although the 
maneuver-induced accelerations were not numerically filtered from the F-106B data, 
they were minimized by data selection. The comparison of normal acceleration data 
is shown in figure 14 and indicates that on the basis of counts per hour the F-106B 
experienced about 10^ more frequently a gust loading of ±0.75g than did the air 
transports. 

A comparison of the derived gust velocities of the F-106B and the air trans- 
ports is shown in figure 15. The L-1011 encountered an updraft producing a derived 
gust velocity almost as large as that encountered by the F-106B (45 ft/sec) but 
downdrafts of corresponding strengths were not encountered by the L-1011. The 727 
encountered downdrafts almost as strong as the F-106B encountered in thunderstorms 
and updrafts 75 percent as strong. The slope of the L-1011 curve with the updrafts 
at high U de values is noted to be about the same as that encountered in the 
thunderstorm flights. The air transports encounter these highest U de only a 
small proportion of the time, but the comparison does show that airliners do occa- 
sionally encounter very turbulent conditions of thunderstorm magnitude. The L-1011 
high turbulence encounter was in the altitude block 9500 to 14 500 ft, which corre- 
sponds to an approach to or departure from the airport where the flexibility to 
deviate around bad weather decreases. The 727 encounter was in the altitude block 
24 000 to 29 500 ft. 
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CONCLUDING REMARKS 


As part of the NASA Storm Hazards Program, loads, motion, and control input 
data were obtained for an F-106B airplane flying inside thunderstorms. The storm 
penetrations were limited to areas of thunderstorms with a radar return of less than 
50 dBZ. Penetrations were made over a pressure altitude range of 10 000 ft to 
38 000 ft with most of the penetrations around 25 000 ft. The data were presented 
in terms of the number of times incremental levels were crossed per hour. The data 
were generally symmetrical about the nominal value. Three different pilots flew the 
airplane for the storm flights and the piloting techniques were basically the same 
for all three pilots. Measured normal acceleration zero crossing rate agrees very 
well with theoretical methods. The maximum derived gust velocity U de of ±45 ft/sec 
measured in this program was about two thirds of the 65 ft/sec measured in the 
Thunderstorm Project. 

The F-106B response to thunderstorms was found to be far more pronounced in the 
roll axis than in pitch. As shown by these data a roll attitude excursion of 10° is 
almost three orders of magnitude more likely to occur than a 10° pitch attitude 
excursion. 

The data were compared with recent airliner data. The comparison showed the 
L-1011 and the 727 had occasionally encountered updrafts or downdrafts of strengths 
almost equal to those encountered by the F-106B (±45 ft/sec) in areas of thunder- 
storms having radar returns less than 50 dBZ. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
March 18, 1986 
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Flight no 


Date 


Pilot 





82-010 

5-23-82 

W. Neely 

82-013 

5-28-82 

W. Neely 

82-017 

6-5-82 

P . Brown 

82-020 

6-27-82 

P . Brown 

82-021 

6-21-82 

W. Neely 

82-023 

6-29-82 

P . Brown 

82-024 

7-4-82 

P . Brown 

82-027 

7-11-82 

P . Brown 

82-028 

7-11-82 

P . Brown 

82-029 

7-16-82 

P. Deal 

82-030 

7-17-82 

W. Neely 

82-032 

7-28-82 

W. Neely 

82-033 

7-28-82 

W. Neely 

82-034 

7-30-82 

P. Deal 

82-035 

7-31-82 

P. Deal 

82-037 

8-6-82 

W. Neely 

82-038 

8-6-82 

W. Neely 

82-039 

8-9-82 

W. Neely 

82-040 

8-9-82 

W. Neely 

82-041 

8-11-82 

P . Brown 

82-042 

8-11-82 

W. Neely 

82-043 

8-17-82 

P. Deal 

82-044 

8-25-82 

W. Neely 

82-047 

9-20-82 

W. Neely 

82-048 

9-20-82 

W. Neely 


Observer 


Location 


P. Deal 

Blackstone, Va. 

B. Fisher 

Annapolis, Md. 

R. Winebarger 

Oregon Inlet, N.c. 

R. Winebarger 

Moyock, N.C. 

R. Winebarger 

Cape Hatteras, N.c. 

B. Fisher 

Nags Head, N.C. 

R. Winebarger 

Lynchburg , Va . 

R. Winebarger 

Blackstone, Va. 

R. Winebarger 

Appomattox, Va. 

R. Winebarger 

Salisbury, Md. 

B. Fisher 

Chesapeake, Va. 

B. Fisher 

Cape May , N . J . 

B. Fisher 

Tangier, Va. 

R. Winebarger 

Easton, Md. 

B. Fisher 

Oceana, Va. 

R. Winebarger 

Smith Point, Va. 

B. Fisher 

Currituck Sound, N.C. 

R. Winebarger 

West Point, Va. 

R. Winebarger 

Lawrenceville, Va. , to 


Franklin, Va. 

R. Winebarger 

Hopewell, Va. 

R. Winebarger 

Scotland Neck, N.C. 

R. Winebarger 

Waverly, Va. 

R. Winebarger 

Warning Area 72A 

R. Winebarger 

Warning Area 72A 

R. Winebarger 

Rocky Mount, N.C., to 


Elizabeth City, N.C. 


















Number of penetrations 








Counts/hr 


t 


All pilots (398.1 min) 



Normal acceleration, g units 

(a) Level crossing rates for all data. 

Figure 4.- Normal accelerations experienced by the F-106B 
while flying inside thunderstorms. 
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All pilots (398.1 min) 
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(c) Peak positive and negative occurrences per storm 
penetration for all data. 


Figure 4.- Continued. 





















Counts/hr 


All pilots (398. 1 min) 
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Lateral acceleration, g units 


(a) Level crossing rates for all data. 

Figure 5.- Lateral accelerations experienced by the F-106B 
while flying inside thunderstorms. 
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Percent of penetrations 



O Pilot 1 (145.4 min) 
□ Pilot 2 (205.4 min) 
O Pilot 3 (47.4 min) 



Lateral acceleration, g units 


(d) Peak positive and negative occurrences per storm 
penetration for each pilot. 

Figure 5.- Concluded. 
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Lateral stick position, percent 


(a) Level crossing rates for all data. 

Figure 8.- Lateral stick movement during thunderstorm penetrations 





Counts/hr 


O Pilot 1 (145.4 min) 
□ Pilot 2 (205.4 min) 



35 39 43 47 51 55 59 63 67 71 75 

Lateral stick position, percent 


(b) Level crossing rates for each pilot. 
Figure 8.- Continued. 


30 



































Roll attitude, deg 


(a) Level crossing rates for all data. 

Figure 9.- F-106B roll attitude during thunderstorm penetrations 
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(d) Peak positive and negative occurrences per storm 
penetration for each pilot. 

Figure 10.- Concluded. 
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O Pilot 3 (47.4 min) 
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Figure 13.- Comparisons of frequency of turbulence encounters 

between three programs . 
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Normal acceleration, g units 


Figure 14.- Comparison of normal acceleration level crossing rates 
for the F-106B while flying inside thunderstorms with three air 
transports in normal service. 
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